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Abstract 

We derive model-independent bounds on production and detection non-standard neutrino inter- 
actions (NSI). We find that the constraints for NSI parameters are around 0(1CP 2 ) to ©(lCP 1 ). 
Furthermore, we review and update the constraints on matter NSI. We conclude that the bounds 
on production and detection NSI are generally one order of magnitude stronger than their matter 
counterparts. 
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I. INTRODUCTION 



The formalism of non-standard neutrino interactions (NSI) is a very widespread and 
convenient way of parametrising the effects of new physics in neutrino oscillations 

y, y, y, Q, y, [9]. Even though present data constrain NSI to be a subleading effect in 
neutrino oscillation experiments, the possibility of their eventual detection or interference 



with neutrino oscillations at present 1(1 



20 



21 



22 



23 



26 



27 



m 12a 



30 



3, Q, Q, H, 



16 



17J and future 



3, 



3l| experiments has triggered a considerable 



interest in the community. In particular, it is a common practice to study NSI in matter, 
which correspond to neutral-current-like operators, assuming that the constraints on the 
NSI affecting production and detection processes are much stronger. However, up to now, 
only model-dependent bounds on such interactions are present in the literature 
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The main aim of this paper is filling this gap and providing model-independent bounds on 
NSI affecting neutrino production and detection processes. Also, the present constraints on 
matter NSI will be reviewed and updated. 

Before entering into details and deriving the current bounds, a discussion on the nat- 
uralness of large NSI is in order. In particular, this argument has been faced for matter 



NSI 



32 



33 



341 ]. but the main message can be applied to production and detection NSI as 



well. Matter NSI are defined through the following addition to the Lagrangian density: 

= -2^f2G F e%\f^Pj\ [V^PlM, (1) 



-NSI 

where / = e,u,d and ejfl encodes the deviation from standard interactions. For example, an 
operator of this kind is induced in fermionic seesaw models once the heavy fermions (singlets 
or triplets) are integrated out leading to a d — 6 operator that modifies the neutrino kinetic 

1. After a transformation to obtain canonical kinetic terms, modified 



energy [35 



36 



37 



couplings of the leptons to the gauge bosons, characterized by deviations from unitarity of 
the leptonic mixing matrix, are induced. Upon integrating out the gauge bosons with their 
modified couplings, NSI operators are therefore obtained. Because of the str ong bounds 
on the unitarity of this matrix, these NSI are constrained to be < O(10~ 3 ) 32j, |37J, |39 |. 



This means that their eventual detection is challenging, although not impossible, at future 



facilities 



40 



41 



42 



43 



44j 



On the other hand, large NSI could be generated by some other new physics, not necessar- 
ily related to neutrino masses, at an energy above the electroweak scale. As a consequence, 



an SU(2) gauge invariant formulation of NSI is mandatory. The simplest gauge invariant 
realization of the operator in Eq. ([1]) implies to promote the neutrino fields to full lepton 
doublets. However, in that case, strong bounds stemming from four-charged-fermion pro- 



cesses would apply 



45|, 



47|. In order to avoid these constraints, cancellations among 



different higher-dimensional operators are required 33|, |45[ . In the case of d = 6 operators 



there is only one combination w 



are also severely constrained 



rich satisfies these conditions and the corresponding NSI 



32 



481 ] . In the case of d — 8 operators it has been shown 



that, avoiding cancellations between diagrams involving different messenger fields or the in- 
troduction of new leptonic doublets that could dangerously affect the electroweak precision 

tests, the only possibilities of evading the constraints imposed by gauge invariance reduce 

I 

to the cases already mentioned, with the consequent stringent bounds [32]. Therefore, in 
order to realise the cancellations that would allow large NSI, some fine-tuning is needed. An 
example of the naturalness prize required is presented in the toy model proposed in Ref. 33]. 
However, even if large NSI are generated in this way at tree-level, dangerous quadratic di- 
vergences contributing to four-charged-fermion operators appear at one- loop 34j. In order 
to have large NSI, another fine-tuning would then be required at one-loop unless the scale of 
new physics is smaller than Attv , where v is the Higgs vev. Alternatively, a symmetry could 
guarantee the cancellation both at tree- and loop-level, but so far no model has been found 
with these characteristics, i.e., leading to large NSI. 

From the previous discussion it is clear that it is not easy to induce large neutrino NSI in 
a specific theoretical framework. However, since it is impossible to exclude them completely 
in a model-independent way and since their effects may be visible at future experiments, we 
think it is worthwhile to derive their present bounds. 

The rest of the paper is organized as follows: First, in Section [Til we define charged- 
current-like NSI and derive bounds on various combinations of e, specifying which bounds 
can be set if only one non-zero e is considered at a time. We then proceed by discussing 
loop bounds on charged-current-like NSI in Section IIHI Finally, we review and update the 
bounds on matter NSI in Section IIVI and make a summary of the results and conclude in 
Section [V] 
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II. CHARGED-CURRENT-LIKE NON-STANDARD INTERACTIONS 



Let us start by considering NSI for source and detector processes. Since these are always 
based on charged-current processes so as to tag the neutrino flavour through the flavour of 
the associated charged lepton, we will refer to these as charged-current-like NSI. The general 
leptonic NSI are given by the effective Lagrange density 

Ci SI = -2V2G F e^ P [! a YPip][^l,PLiys], (2) 

where P is either Pl or Pr and, due to Hermiticity, e"g P = e^ P * '. For charged-current-like 
NSI a 7^ (3; in particular, a = fi and (3 = e are the only parameters of importance for 
neutrino oscillation experiments due to their effect in neutrino production via muon decay. 
Notice that a = (3 = e would instead correspond to matter NSI. 

In a similar fashion, the charged-current-like NSI with quarks are given by the effective 
Lagrange density 

C q mi = -2V2G F e% P V qql [qYPq'}[Ll,P L ^} + h.c, (3) 

where q is an up-type and q' is a down-type quark. Naturally, only q = u and q' = d are 
of practical interest for neutrino oscillations, due to their contributions to charged-current 
interactions with pions and nuclei. Because of this, we will concentrate on constraining e^a 
as well as e^l. Since the relevant combinations of NSI that contribute to some processes will 
be of an axial or vector structure we define 

-y&V _ 7<5_R , 7<5L ( A \ 
t a/3 ~ t a(3 ' t af3 ' 

~{8A _ ~{8R _ -ySL /r\ 
t a/3 ~ t a(3 t a/3 ' \° J 

in order to simplify the notation. Notice that more general Dirac structures such as scalar 
or tensor couplings can in principle be considered to generalise Eqs. ()2]) and ([3]). However, 
these NSI will have the wrong chirality to contribute coherently with the SM production 
and detection processes -something that is usually assumed for NSI- and therefore linear 



interference of these NSI will require an extra chirality suppression 25|, |49(. For this reason 
we will neglect them here. 
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A. Bounds from kinematic Fermi constant 



At present, the most precise determination of the Fermi constant Gf is through the muon 
decay rate. However, if NSI of the form e^g are present, this will make the measured Fermi 
constant from muon decays differ from the true Fermi constant according to the relation 
Gfj, = Gpf(e^ L , ^ Q) 3p k«/3 P | 2 )- Here we have introduced the function 

f(x,y) = l + 2Re(x)+y, (6) 

where x represents the interference between the SM and the particular NSI that contributes 
coherently with the SM to the process and y is the incoherent sum of the NSI contributions. 
Therefore, in all the processes considered, stronger bounds will be implied for the real part 
of x. Given the relation between G^ and Gp, an independent measurement of the Fermi 
constant will constrain e^jf . We will consider two different ways of deriving the value of Gf, 
one involving only the kinematic measurements of the gauge boson masses and one involving 
comparison to the quark sector. 

For determining Gf from kinematic considerations, we need to review the predictions of 



the Standard Model. From Ref. 



501 ]. we have 



M W = (7) 

where A = \J Tca/(y/2Gp), = 1 — M^/M§, a is the fine-structure constant, and Ar = 
0.03690 ± 0.0007 is the radiative correction to the tree-level relation. Thus, we obtain the 
relation 

G %aM2z (8) 

F y/2M* r {M*-M* r )(l-Ary 

For the masses of the vector bosons, we use the combined fit for the W mass from LEP and 



50]. 



Tevatron, M w = 80.398 ±0.025 GeV, as well as M z = 91.1876 ±0.0021 GeV from LEP 
The resulting Fermi coupling constant is 

Gf = (1.1696 ± 0.0020) • 10~ 5 GeV~ 2 (la). (9) 

Comparing with G M , we obtain 

% = m; L , E = 1 16637 ±a ° 0001 = 0.9973 ± 0.0017, ( io) 

G F " ^ 010 1.1696 ±0.0020 V ; 
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which represents a 90 % confidence level agreement with the Standard Model expectation. 
The only truly model-independent bound that we can extract from this is on the combination 
f{ £ e^, L y S a /3P | 2 )- O n the other hand, it is common practice to assume the presence of 
only one non-zero £ at a time in order to avoid cancellations inside f(£^f',^2 a pp \ £l £p > \ 2 )- 
In this way, the following bounds can be obtained: 1 

Re«0 = (-1.4±1.4)-10- 3 , (11) 
IC/I < 0.030, (12) 

at 90 % confidence level. 



B. Bounds from CKM unitarity 

One way of constraining the completely leptonic NSI, as well as some of the charged- 
current NSI with quarks, is to make the assumption that the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix is unitary, as predicted by the Standard Model. The experimental test of the 
CKM unitarity is essentially based upon the determination of V u d and V us from beta- and 
Kaon-decays, 2 where the Fermi constant extracted from muon decay is used to predict 
the decay rates. These are proportional to 

T ex G 2 F \V UX \ 2 , (13) 

which means that, by inserting in place of Gf, we are actually determining |V^f| 2 = 
a/3P \ £ a/f\ 2 )- Adding the information from beta- and Kaon-decay experi- 
ments and assuming that leptonic NSI dominate over quark NSI, we have 50j 



\V5 1 2 + |K*T = \YuA±SY^t = \ = 0.9999 ± 0.0010 (14) 

J l fc eM J l^aQP l fc a/3 I ) J l fce M ' Z-^a(3P l fc a/3 



at lcr, where the CKM unitarity is inserted in the second step. Again, this translates into a 
bound for f(s^ L , J2 a /3P ka/Tl 2 ); ^ut making the assumption of having only one non-zero e 



at a time we obtain: 



Re(0| < 4.0 -10- 4 , (15) 



1 Throughout the paper we will follow the statistical approach proposed in Ref. [5l| by Feldman and Cousins. 

2 In principle V u b should also be considered. However, its value is smaller than the uncertainty in the other 
two matrix elements and we therefore leave it out of our discussion. 
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IC/I < 0.030, (16) 

at the 90 % confidence level. Notice that the bound of Eq. (|T5|) is slightly stronger than the 
one obtained from the kinematic determination of the Fermi constant, but it relies on one 
extra assumption, i.e., the unitarity of the CKM matrix. 

On the other hand, if we assume that the NSI with quarks are dominating, then the 
insertion of G M in place of Gf is not leading to any ambiguities. However, NSI of the form 
e ud will contribute to the beta-decay rate, through which V ud is extracted. Experimentally, 
only superallowed + — > + decays are considered, which means that the nuclear matrix 
element will have a vector structure and, therefore, only the vector NSI combination will 
contribute in the following way: 

r^G F \V ud \ 2 f(e u J v ,J2\^ V \ 2 )- (17) 

a 

Since the Kaon decays are not affected by e ud , these can be used to extract V ud indi- 
rectly from the assumption of CKM unitarity {i.e., \V u d\ 2 = 1 — |V„ S | 2 ). The result of 



this operation is \V ud \ 2 = 0.94915 ± 0.00086 50), which should be compared to the value of 



\V ud \ 2 = \V ud \ 2 f{e^ v , J2 a ke^H derived from beta decays \V ud \ 2 = 0.94903 ± 0.00055 |50| 



Once again a truly model-independent bound can only be extracted for the combination 
f( £ ee V y ^ Q l £ ea V | 2 )' t» u t making the assumption of taking one e at a time we obtain: 

\Re(e^ v )\ < 0.00086, (18) 
\e^ v \ < 0.041. (19) 

Notice that, unlike the determination through the kinematic Gp, the determination of 
the non-standard parameters ef^L through CKM unitarity relies on the assumption that 
the quark interactions are not affected, making the resulting bounds slightly more model- 
dependent. On the other hand, if a given model predicts both lepton (s^V) and quark 
(e^l) NSI simultaneously, the bounds on from the kinematic Gf compared to muon 
decay would still apply, while somewhat weaker bounds on e™| could still be derived after 
propagating the errors derived on the former through the CKM unitarity relation. 
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C. Bounds from pion processes 



For the quark charged-current NSI involving charged leptons other than electrons, the 
universality tests stemming from the relative decay rates of charged pions as well as that 
of taus into pions can be used to set bounds. The squared and summed matrix element 
involving a charged pion, a charged lepton and a neutrino is modified according to 

£ \M{*, 4, ^)l 2 = \m(tt, e a , v a )\ 2 fKi A , £ KT\ 2 )- (2°) 

p & 

This modification is equivalent to violations of weak interaction flavor universality identifying 
9 2 f{ eU aai Xw? ka/? 4 ! 2 ) = 9ai where g a is the W coupling to the lepton flavour a. Comparing 
the rates of tt — > eu, tt — > and r — > ttu, bounds can be set on the ratios g a /gp. From 
Ref. (52! we have 

^ = 1.0021 ±0.0016 and — = 1.0030 ± 0.0034 (21) 

9e 9 pi 

at la. Thus, if only one e is considered at a time, we obtain the following bounds at the 
90 % confidence level: 

(22) 
(23) 
(24) 
(25) 
(26) 
(27) 

Notice that the bounds on \e^ A \ are more stringent than the bounds on because the 

offset of the best-fit from the Standard Model expectation goes in the opposite direction 
with respect to the effect of 

It is important to note that a model that predicts equal f{s u ^ A ^ J2p \ £ a$ A \ 2 ) f° r a = e, /1, r 
cannot be bounded using this type of argument, since it affects all of these decays in the 
same way and universality is not violated. However, if we only consider e of one chirality 
at a time, then this would imply that e^ p and e^ p share the stronger bounds derived for 
e^f from the CKM unitarity. 



Re(e^) 




(2.1 ±2.6) ■ 10~ 3 , 


I p udA 1 


< 


0.078, 


Re{e u T dA ) 




(3.0 ±5.5) • 10" 3 , 


I p udA 1 


< 


0.13, 


Re(e^) 




(-2.1 ± 2.6) ■ 10~ 3 


1 udA 1 


< 


0.045. 



s 



In a similar fashion, we can use the universality test between the /i — > evv and r — > 
yuz/z/ decays to constrain the non-standard couplings e^. This constraint is related to the 
lepton universality ratio g T /g e = 1.0004 ± 0.0022. Therefore, the inverse of this number 
is a measurement of f(eefi L , Y2 a /3P \ e ^p\ 2 )i where we disregard possible modifications of 
the tau decay which are not important for neutrino oscillation experiments. The resulting 
bounds are: 

Re(^) = (-0.4 ± 3.5) • 10" 3 , (28) 
\e^\ < 0.080. (29) 



D. Bounds from oscillation experiments 

Production and detection NSI imply that a neutrino produced or detected in association 
with a charged lepton will not necessarily share its flavour. This means that flavour conver- 
sion is present already at the interaction level and "oscillations" can occur at zero distance. 
Indeed, in the presence of NSI, 

P aP {L = 0) ~ |e#f (30) 
if the neutrino is produced through pion decays and 

P ea (L = 0) ~ £ \eff\ 2 as well as P M/3 (L = 0) ~ £ \e»f\ 2 (31) 

j3P aP 

for neutrinos produced through muon decays. For the detection through inverse beta decays 
the situation is a bit more involved since the relative contributions of the different chiralities 
vary depending on the energy regime due to the nuclear matrix elements. Here we will 
discuss the cases of very low (E < 1 GeV) and very high (E > 10 GeV) energies. In the 
first case the neutrino-nucleon cross section is proportional to (g v + 3g\), where gy — 1 and 
gA = 1.23. This means that the vector and axial combinations of the NSI that can mediate 
the processes will contribute incoherently with those relative strengths to give: 

Pafi(L = 0) ~ Y^jd^T + 3fli|e£T). (32) 
Notice that, if only one non-zero e with definite chirality is present, then 

P a ,(L = 0) ~ \eff | 2 . (33) 
9 



Experiment 


Channel 


Bounds 


KARMEN 


Vfi v e 


iieP 

toe 


< 0.025, 


p udA 
^efi 


< 0.028, 


udV 
'-■efi 


< 0.059 


NOMAD 


-> U T 


1 p udA 


< 0.013, 




< 0.013, 


p udR 

t-Tfl 


< 0.018 


NOMAD 


v e -> V T 


aeP 

C-CtT 


< 0.087, 


1 fc re 


< 0.087, 


\ p udR 
| fc Te 


< 0.12 


NOMAD 


Vfi -> v e 


1 F udA 
\ fie 


< 0.026, 




< 0.026, 


pudR 


< 0.037 



TABLE I: Bounds (90 % CL) from oscillations at zero distance. In each line, the first bound refers 
to production NSI, while the other two are for detection NSI. 



We will make this assumption when we will summarise the bounds in the last Section. On 
the other hand, at very high energies, in the deep inelastic scattering regime, the left-handed 
NSI contribute to the neutrino cross-sections with a strength about twice that of the right- 
handed, the actual factor being given by the ratio of the neutrino and antineutrino cross 
sections at high energies for an isoscalar target r = a u /a D ~ 6.7/3.4 = 1.97. We then obtain: 

P a p(L = 0) ?± \ef a L | 2 + ~K*| 2 . (34) 

We can therefore use the very precise constraints on flavour oscillations from experiments 
such as KARMEN 53J and NOMAD 5J, |55j. Motivated by the large mass hierarchies 
and small mixing angles observed in the quark sector, these experiments explored neutrino 
oscillations at very short baselines with high precision and no evidence of flavour change 
was found. Both KARMEN and NOMAD produced neutrino beams from 7r + decays as well 
as the subsequent fi + decays and detected them through inverse beta decay. In the case of 
KARMEN the neutrinos were produced via /i decays at rest, so that the neutrino energy was 
always below 50 MeV. On the other hand, NOMAD aimed at the detection of z/ T , so higher 
energies ~ 20 GeV were exploited. Table [T] contains a summary of the different oscillation 
channels they explored and the bounds they imply for the NSI parameters. 



III. LOOP BOUNDS 



The tree level effects of neutrino NSI are difficult to constrain since neutrino detection 
and flavour tagging is challenging. However, NSI may mix with four-charged-fermion op- 
erators at the loop level inducing flavour-changing charged-lepton interactions, for which 
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FIG. 1: (a) The vanishing one-loop contribution to the mixing in the running between the matter 
NSI and the four-charged-fermion operator via W exchange, (b) The non-vanishing one-loop 
contribution to the mixing in the running between the charged-current-like NSI and the operator 
inducing /j, — > e conversion in nuclei. 



strong bounds exist. In Ref. 34j] it was shown that, for a certain class of diagrams (see 
Fig. [T^,), the logarithmic divergences that would indicate the mixing in the running between 
NSI and four-charged-fermion operators canceled. Therefore, only model-dependent finite 
contributions remain and no model-independent bound can be derived through one-loop 
considerations. We have checked that this is also the case for most neutrino NSI at produc- 
tion and detection. There is, however, an exception: the NSI parameter mixes with the 
operator that induces muon to electron conversion in nuclei through the digram of Fig. [TJd. 
The computation of this diagram yields a logarithmic divergence: 

A 



3V2G F ael d e L , 

-log 



[u~fpP L u][nYP L e]. 



(35) 



27TS* ~~~\M Wi 

Since the coefficient of this divergence can be interpreted as the coefficient of the logarithmic 



running of this operator, we can estimate the bound by assuming that log(A/M| 



gives a contribution to /i — > e conversion in nuclei of the form (see, e.g., Ref. 561]): 



w 



R(/x~ -> e~) 



m 5 J2VW + V^f\C\ 



T(n capture) 

where C is the coefficient of the operator in Eq. (1351) . Using R(/x 



1. This 



(36) 



< 7.0 • 10" 13 for 



conversion in Au 



50j as well as = 0.0974 and V^™ = 0.146 |56j, a very strong bound 
on the NSI is derived: 



\e u ^ d e L \ < 1.8 • 10" 



(37) 



We would like to remark that, also in this case, a quadratic divergence is present. In principle, 
this contribution could dominate over the logarithmic one, but its value is model-dependent 
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and reliable bounds cannot be derived from it. The contribution from the logarithmic run- 
ning could also be canceled, but only at a given scale, which makes the resulting constraint 
more reliable. 



IV. NEUTRAL-CURRENT-LIKE NON-STANDARD INTERACTIONS 



For completeness, we will now also review the current status of the bounds on NSI 
matter effects, or neutral-current-like NSI defined in Eq. (OQ). This type of NSI is the most 
extensively studied in the literature, since it has been generally assumed that the constraints 
on the charged-current-like NSI are much stronger. We would like to stress that, in specific 
models, charged-current-like and neutral-current-like processes are expected with similar 



strengths 



32] 



In most phenomenological studies the NSI parameters are reduced to the effective pa- 
rameters 



E 



(38) 



where rif is the number density of the fermion /. This is the natural parameter in neutrino 
oscillation analyses since it corresponds to the replacement 



matter 



V" 



V 



1 





\ 



V 



/ 



\ 



1 





\ 



^■ee ^-efi ^er 



\ 



e e^i £ W 



J] 



(39) 



in the matter interaction part of the neutrino flavour evolution. Thus, assuming uncorrelated 
errors, the bounds on e a/ 3 could be approximated by 

1/2 



^< E[K) 2 + ( 3 ^) 2 + ( 3 « ; 



(40) 



for neutral Earth-like matter with an equal number of neutrons and protons and by 

1/2 



(41) 



for neutral solar-lf 



from Refs. [53, 



ce matter, consisting mostly of electrons and protons. Using the bounds 



58 



59 



601 ] . but discarding the loop constraints on [34j . the resulting 
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bounds on the effective NSI parameters would be 

/ An nQQ q n \ 

and 



l fc a/3l < ~ 



4.2 0.33 3.0 
0.33 0.068 0.33 
3.0 0.33 21 



r° I < 



2.5 0.21 1.7 
0.21 0.046 0.21 
1.7 0.21 9.0 



(42) 



respectively. Notice that atmospheric neutrino oscillations also constrain the values of matter 
NSI through the relation e\ 



. 7T ~[|^ 2 ±0(O/L)]/(l + ej 
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significantly smaller than one, this would set a stronger bound e® T < O(10) 



6l| . As long as 1 + £® is not 

< 

We want to stress the fact that the constraints on e e , e u and e d have been derived under 
the assumption of taking one non-zero e at a time. Thus, the approach of combining them 
together as in Eq. (j42p is not fully consistent. For this reason, in the compilation of all the 
results in the following Section, the bounds will be quoted separately. 



V. SUMMARY OF RESULTS AND CONCLUSIONS 



In order to easily overview our results, we here present the constraints from the previ- 
ous sections in tabularized format. In Tab. [IT] we present the different bounds available for 
£ a3 wnne the bounds for are presented in Tab. IIHI Taken all together, the most strin- 
gent bounds available for both charged-current-like and neutral-current-like NSI relevant for 
terrestrial experiments are given by: 

/ n noc; n dqo n nm \ 



-.ud I 

-a/3 1 



-a/3 1 



< 



< 



< 



0.025 0.030 0.030 
0.025 0.030 0.030 
0.025 0.030 0.030 



/ 



0.041 


0.025 


1.8 • 10~ 6 
0.026 


0.078 


0.087 


0.013 


0.12 


0.018 


0.06 1Q 
0.14 U U 


0.4 
0.27 


0.10 0.03 


0.10 


°- 4 10 
0.27 


0.16 
0.4 



0.13 



(43) 



(44) 



(45) 
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ue 


Kin. G F (L,i?) 


CKM unit. (V) 


Lept. univ. (A) 


Oscillation (L, R) 


ue 
&ee 


< 0.030 


< 0.030 


< 0.080 


< 0.025 


ue 
&eu 


(-1.4 ± 1.4) • 10" 3 (]R,L) 
< 0.030 


< 4- 10~ 4 (M) 
< 0.030 


(-0.4 ±3.5) • 10~ 3 (R) 
< 0.080 




ue 
ter 


< 0.030 


< 0.030 


< 0.080 


< 0.087 


ue 
^{xe 


< 0.030 


< 0.030 


< 0.080 


< 0.025 


ne 
MM 


< 0.030 


< 0.030 


< 0.080 




lie 

t-LIT 


< 0.030 


< 0.030 


< 0.080 


< 0.087 


ue 


< 0.030 


< 0.030 


< 0.080 


< 0.025 


Lie 


< 0.030 


< 0.030 


< 0.080 




ue 


< 0.030 


< 0.030 


< 0.080 


< 0.087 



TABLE II: Bounds (90 % CL) on the purely leptonic charged-current-like NSI e^, relevant to the 
neutrino production through muon decay, e.g., at a Neutrino Factory. The letters L,R,V,A refer 
to the chirality of the e which is actually bounded, while R stands for the real part of the element 
only. See the text for details. 



/ 1.0 
0.7 



\ £ af3\ < 



0.05 0.5 



0.05 0-003 0.05 



0.5 0.05 



\ £ al3\ 



< 



V 

( 0.3 
0.6 

0.05 



(46) 



1.4 
3 / 
\ 



0.05 0.5 
0.05 



0.003 
0.015 



V 



0.5 0.05 y 
6 



(47) 



/ 



Here, whenever two values are quoted, the upper value refers to left-handed NSI and the 
lower to right-handed NSI. We would like to stress that, before applying these constraints, 
the reader should refer to the appropriate Sections in order to be aware of the assumptions 
under which they were obtained. 

To summarise, we have presented the model-independent bounds that can be derived 
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-.ud 

£ a/3 


CKM unit. (V) 


Lept. univ. (4) 


Oscillation 


Loop (L) 


pud 


< 8.6 • 10~ 4 (R) 
< 0.041 


(-2.1 ± 2.6) • 10~ 3 (R) 
< 0.045 






p ud 


< 0.041 


< 0.045 


< 0.028(4) 

< 0.059(V) 

< 0.026(L) 

< 0.037(1?) 




gud 


< 0.041 


< 0.045 






^ud 




< 0.078 


< 0.026(4) 


< 1.8- 10~ 6 


p ud 
MM 


- 


(2.1 ±2.6) • 10~ 3 (R) 
< 0.078 


- 


- 


p ud 

fc /IT 


- 


< 0.078 


< 0.013(4) 


- 


piid 

TG 




< 0.13 


< 0.087(L) 

< 0.12(12) 


- 


£.ud 
r M 




< 0.13 


< 0.013(L) 

< 0.018(12) 




p ud 

TT 




(3.0 ±5.5) • 10~ 3 (R) 
< 0.13 







TABLE III: Bounds (90 % CL) on the quark charged-current-like NSI e^, relevant to the neutrino 
production through hadron decays as well as detection processes. The letters L, R, V, A refer to 
the chirality of the e which is actually bounded, while R stands for the real part of the element 
only. See the text for details. 



for various types of NSI. Since the neutral-current-like NSI have been studied extensively 
in the literature and the bounds on these are fairly well known, we have just summarised 
these results and concentrated on the charged-current-like NSI, which usually are simply 
considered to be very strongly bounded, although no model-independent analysis has been 
readily available. The result of our analysis is that the charged-current-like NSI, which 
are of interest mostly for their impact on neutrino production and detection, are generally 
bounded by numbers of (9(10~ 2 )-(9(10 _1 ), except for the very strong loop bound on e u ^ L 
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due to the operator mixing inducing /j, — > e conversion in nuclei. We find that these bounds 
are about one order of magnitude stronger than the bounds on the neutral-current-like NSI. 
We therefore argue that production and detection NSI should not be neglected with respect 
to matter NSI, especially taking into account that, in most realisations, both kinds of NSI 
are induced with similar strengths. Moreover, NSI saturating the bounds derived here will 
be within the sensitivity reach of planned neutrino oscillation experiments. However, as 
discussed in the introduction, most models leading to NSI generally affect other processes 
and therefore stronger bounds than the ones derived here apply. 
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